Are cloned quantum states macroscopic? 
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We study quantum states produced by optimal phase covariant quantum doners. We argue that 
cloned quantum superpositions are not macroscopic superpositions in the spirit of Schrodinger's cat, 
despite their large particle number. This is indicated by calculating several measures for macroscopic 
superpositions from the literature, as well as by investigating the distinguishability of the two 
superposed cloned states. The latter rapidly diminishes when considering imperfect detectors or 
noisy states, and does not increase with the system size. In contrast, we find that cloned quantum 
states themselves are macroscopic, in the sense of both proposed measures and their usefulness in 
quantum metrology with an optimal scaling in system size. We investigate the applicability of cloned 
states for parameter estimation in the presence of different kinds of noise. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.65.Yz, 03.65.Ta 
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Quantum mechanics predicts puzzling features of 
quantum systems, mainly due to the possibility of su- 
perposition and interference. While such properties are 
regularly observed at a microscopic scale, the produc- 
tion of large-scale entanglement or superpositions of 
macroscopic objects has become the focus of interest 
in recent years. Remarkable experimental progress has 
been reported with various setups, including atomic sys- 
tems [TJ , superconducting devices [2] [5] , and optome- 
chanical systems |3] as well as photons |BJ ■ Some of 
these experiments aim for the production of a quantum 
superposition of macroscopic objects, in the spirit of the 
famous Schrodinger cat gedankenexperiment |7]. 

However, the required features that make a quantum 
object macroscopic or a macroscopic superposition are 
not obvious. Should one call a single atom in a su- 
perposition of two localized states macroscopic if the 
two wavepackets are far apart? What about massive 
objects, or objects consisting of large number of parti- 
cles, where the displacement of the two wavepackets is 
extremely small? Several attempts have been put for- 
ward in the literature that aim for providing measures of 
macroscopicity and macroscopic superpositions |8"rtT5], 
Here we will concentrate on discrete systems consist- 
ing of N qubits, where the notion of macroscopicity is 
better explored [T5]. In this case it is evident that a 
large particle number N is required to call such a state 
macroscopic; however, this is apparently not sufficient. 
For instance, a product state of N qubits is certainly not 
macroscopically entangled, and also a quantum state 
where two particles are in a superposition but the rest 
is in some fixed product state is not a macroscopic su- 
perposition either. 

We consider quantum states produced by an optimal 
phase covariant quantum doner (PCQC) [IB]. These 
states have been considered in the context of photons, 
where it has been demonstrated that an optimal PCQC 
can be experimentally realized by using parametric am- 
plification [17], leading to quantum states consisting of 
tens of thousand photons (5J |BJ HH1 HHI • When applied 
to an initial state in a superposition, or part of a max- 
imally entangled photon pair, this leads to quantum 
superpositions of two such multiphoton states. It has 
been argued that with such states it is possible to ob- 



serve entanglement with the naked eye |18] and that in- 
deed this corresponds to the production of micro-macro 
entanglement in the spirit of Schrodinger's cat [SUB]. 

Here, we critically investigate the macroscopicity of 
such systems. Thereby, we carefully distinguish be- 
tween macroscopicity of the states as such and macro- 
scopic superpositions |15| . By applying different mea- 
sures from the literature |llffT51 [T5] that have been pro- 
posed to determine an effective size of quantum super- 
positions, we show that all these measures render cloned 
superpositions microscopic, despite their large particle 
number. This conclusion is also supported by consider- 
ing the distinguishability of the two superposed states, 
which rapidly diminishes when considering noisy states 
or imperfect detectors (see also Refs. |18'.f20]) and which 
does not increase with the number of particles. In con- 
trast, cloned states as such turn out to be macroscopic, 
both in the sense of other proposed measures [101 [T5] , 
but also with respect to their applicability for parameter 
estimation with enhanced precision, since the Heisen- 
berg limit can be reached [21] . We investigate the ap- 
plicability of cloned states for quantum metrology |55] 
in the presence of different kinds of noise and imperfec- 
tions. 



Phase covariant quantum doner 

We start by reviewing the optimal PCQC [TB]. We 
consider systems consisting of N qubits with Hilbert 
space H = (€ 2 ) <g > N . A PCQC produces from a single 
qubit state \(p) = -t=(|0) + e lv \l)) lying on the equa- 
torial plane of the Bloch sphere N optimal copies, in 
the sense that the resulting symmetric iV-qubit state 
is such that the fidelity of each single particle reduced 
density operator with respect to the initial state is max- 
imal. The optimal unitary cloning operation performs 
for odd N the mapping [TB] 



|±) -► |^) 



V2 



N, 



N -1 



± 



*^±i'>l.(l) 



where |±) = -4= (|0) ± |1}) are the eigenstates of a x , and 
\N,k) denotes a 7V-qubit symmetric Dicke state with 



k excitations [23] , that is, a state given by a coherent 
superposition of all possible permutations of product 
states with k ones and N — k zeros. Note that, for even 
TV, Eq. (Ill) differs in the details. However, the results 
and conclusions are equivalent to odd TV, on which we 
focus for the remainder of this article. 

Applying the cloning transformation to the second 
qubit of an entangled state -L(|+) (g) |— ) — |— ) <g) |+)) 
results in an TV + 1 qubit state 
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This is the state we will consider in the following. Notice 
that it was argued in Refs. [HI H3 [H] that Eq. Q cor- 
responds to the quantum state that is generated in a 
photonic setup using parametric amplification. 



Measures for macroscopic superpositions 

We now review and apply some measures to quantify 
the "effective size" N e s of multipartite spin states of the 
form 
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(I0o) + I<M). 
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For the state (pi , we use the splitting 
and \<j)i) = \-)® |V> + >. 

The first two measures try to capture the notion of 
"macroscopically distinct" states |0q) and \<j)\). Ko- 
rsbakken et al. [H] argue that these two states are 
macroscopically distinct if they are distinct already on 
a microscopic subset of the whole number of particles. 
Hence the effective size TVjg of a quantum state (pi is 
defined as the maximal number of partitions among the 
particles such that measuring just one of these groups 
allows us to differentiate between \(f>o) and \<f>i) with 
high probability. 

To apply the measure we first note that the "mi- 
cropart" in (pi constitutes one group which allows with 
a o~ x measurement a perfect distinction between \<f> ) 
and |</>i). The nontrivial part is to determine the opti- 
mal division of the doner states with respect to a suc- 
cess probability P for their distinction. This proba- 
bility is defined as P := \ + \\p\ - pi ||i [25J with 
Pk = Trjv\ fe |V' ± )(V' ± I- The detailed calculation is pre- 
sented in the Appendix (see also Ref. [1]). For k = 1, 
we have P = |[1 + 0(1/TV)]. If we increase k, the suc- 
cess probability goes up, but saturates to a value below 
0.82; see Fig. [T] On the other side, we can also measure 
all but k particles. From Fig. [l] we see that it makes 
little difference whether we measure 100 or TV— 100 par- 
ticles for TV ^ 1. Even if we measure TV — 1 particles, 
P « |(1+^ts) ~ 0.85. The effective size depends on the 
uncertainty in the measurements we are willing to allow. 
Nevertheless, it seems clear that we have to measure all 
particles to be certain whether we encounter IV* ) or 
|i/)~). Hence the state is not a macroscopic superposi- 
tion; that is, TV o ^(V0 = 2. 
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FIG. 1: Success probability P to distinguish \ip + ) from 
\ip~) by measuring k out of TV particles. The results are 
valid in the limit of large TV. We clearly see that there is 
little difference between the scenarios where one measures 
100 particles and all but 100 particles, for which we have in 
both cases P « 82%. 



Similarly, Marquardt and co-workers [13] see two 
states distinct on a macroscopic scale if one has to ap- 
ply on average O(N) one-particle operators on |</>o) in 
order to achieve a high overlap with \<f>i). For the state 
\ip) of Eq. (pi, we need at least two one-particle oper- 
ations because the microscopic as well as the macro- 
scopic parts are orthogonal on their own. But we will 
see that we also do not need more than two; therefore 
N^g(ip) — 2, which means that the superposition is not 
macroscopic in this sense. We prove the statement with 
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the operator a z ®M z , where Mj = J2 
|26| . The Dicke states are eigenstates of M z with 
M Z |TV, i(TV±l)) = =F |iV, |(TV± 1)), which means 
that a z (g) M z |0o ) == \4>i)- With a z ® M z as a linear 
combination of two-particle operations, we entirely map 
|0o) onto \<fii) and vice versa. 

Bjork and Mana [11_ call a superposition of the type 
(pi macroscopic if the usefulness of \(f>) is largely in- 
creased compared to the single constituents \<J)q) and 
\4>i). This usefulness was defined in terms of interferom- 
etry experiments. For this it suffices that a local Hamil- 
tonian H generates a unitary evolution which is much 
more rapid for \4>) than for \<fr ) and \(j)\). The "speed of 
evolution" is measured by the time 8± it takes to evolve 
the state to an orthogonal one. The effective size of 
the state Q is hence the ratio TV^(0) = 6±(4>o)/6x(4') 
|27| . Since for spin systems TV^j is not always defined 
(e.g., when \<f> ) is an eigenstate of H), we suggest to use 
instead the "relative Fisher information" [TS], which is 
based on the quantum Fisher information J-(p, H) [55] 
(see also the last part of the Appendix) . Note that for 
pure states \<j>), T(4>,H) is proportional to the variance 
of H , that is, J-{(j>, H) = AV^H). To measure the use- 
fulness of \<p) compared to its constituents in Eq. (pi, we 

set TV c r |(<T)) = ms^H-.iocaiV<t,{H)/[mzx Ho docaiV<t>a{Ho)]- 
If N*g(<t>) = O(N), we call \4>) a macroscopic superposi- 



tion due to Bjork and Maria; a constant scaling signifies 
a "microscopic" superposition. We hence calculate the 
maximal variance of \ip) and \4>q) with respect to all 
sums of local operators with a constant operator norm. 
For Dicke states with N/2 excitations this has already 
been done in Ref. [21 , by using the so-called covari- 
ance method. For the micro-macro state \ip) this can 
be easily adapted (for details, see the Appendix). The 
optimal operators yield V^,(ff opt ) = §(iV+l) 2 +iV+l 
and V 0o (#o pt ) = |(JV + l) 2 . The ratio of those two is 
the relative Fisher information N^(ip) = l + 2/(N + l). 
We see that the quantum state is a microscopic super- 
position, because the two constituting states \<J)q) and 
\4>i) are very similar to the total state \ip) with respect 
to the maximal variances. 



a = y/N, but the conclusions are the same for con- 
stant a > 1. Last, we investigate D if the states under 
consideration are subject to noise. We study a local 
phase-noise channel £ z . This map acts on an arbitrary 



N 



input state p as £ z (p) — Yii=i up + (I — u)a z pa z with 
the noise parameter u — i(l + e~ 7 *) £ [^ , 1] . The 
corresponding A Noisc equals EjTr [in£ z (|V> + X^ + I)] ~ 
Tr[7r^ z (|V-)(V-|)]|- 

The numerical calculations are presented in more de- 
tail in the Appendix and, for a certain choice of the 
parameters, are plotted in Fig. [2] For all three scenar- 
ios of imperfect measurements or state preparations, a 
numerical extrapolation shows that the distinguishabil- 
ity between \ip + ) and \ij}~) converges to the same value 
D ss 82% [30J and cannot be improved by increasing N . 



Distinguishability of cloned superpositions 

In Refs. [5J El HI] it has been argued that superposi- 
tion (pi is macroscopic since the two macro states \tp + ) 
and l^ - ) are "macroscopically distinguishable." This 
characterization is footed on the difference of the expec- 
tation values for the magnetization in the x direction, 
which scales linearly with N, since (M x )^± — ±\N . 
This reasoning implicitly assumes that the two states 
l^) are better and better distinguishable as N grows, 
especially if noise or imperfections of the measurements 
are taken into account. 

Here we show that this is not the case (see also 
Refs. [18, 20 ). If we assume noiseless measurements, 
the two states are perfectly distinguishable with respect 
to M x for all N, because l^) live exclusively in sub- 
spaces of M x with an even and odd excitation num- 
ber, respectively. On the other side, the variances of 
IV 7 ) with respect to M x scale quadratically with N, 
Vjij±(M x ) = \ {N 2 — l), which hinders an arbitrarily 
high distinguishability in the presence of noise even in 
the limit of large N. Given the spectral decomposition 
of a measurement M = y\ TOjTTi with the eigenvalues 
mi and the (possibly high-dimensional) projections ttj, 
the quantity D = \ + \& with A = J^K 71 *)^ —(^i)il>- I 
serves as a characterization of the distinguishability be- 
tween \ip + ) and \ip~)- 

Here, we consider three different scenarios for M. 
The first scenario is a toy model which mimics a non- 
perfect resolution of a measurement apparatus. We 
therefore group the N + 1 eigenvalues of M x into 
pairs, such that m,, rrii+i — > | (toj + Jnj+i) if i is odd. 
This means that we cannot distinguish between even 
and odd parity of the M x measurement. Hence, one 

has A Pair = Ei=i,3,...l( 7r i + ^+1)0+ - (tt« + Ti+i>^-|- 
The second approach is to regard a generalized posi- 
tive operator- valued measure (POVM). The measure- 
ment outcomes mi are associated with the Kraus op- 
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erators E± — 127=0 a/"-7 6 ^ ^ 4<T ^-i where n, 



J2k=o e _(fe_i) /(2ct ) is the normalization of the POVM. 
The parameter a > reflects the normally distributed 
uncertainty induced by the measurement. Hence we 
study A POVM = EilC^f )v+ - (£?),/,- |. Here we choose 
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FIG. 2: Distinguishability D for the three scenarios of 
imperfections discussed in the text. The plot indicates that 
increasing the system size does not improve the success prob- 
ability to distinguish \ip + ) from \tp~)- 

We conclude that -according to the literature |1 H — 
\T3\- and arguments based on distinguishability, one 
cannot use optimal covariant cloning devices to gen- 
erate macroscopic superpositions. 



Usefulness for parameter estimation 

In contrast to measures on macroscopic superposi- 
tions |llllT3] there exist also proposals for a charac- 
terization of macroscopic quantum states that do not 
rely on a certain structure of the state as in Eq. pi 
|10l [21 US]- Directly applicable to spin systems are 
|101 ITS] , which qualify the cloned superposition as 
macroscopic. 

The so-called index p |10| -introduced by Shimizu 
and Miyadera- formalizes the observation that classi- 
cal probability distributions encountered in statistical 
mechanics show sharp peaks for local observables M 
like the magnetization. This comes from the fact that 
these M exhibit an expectation value (M) oc N and a 
variance V(M) oc N, which means that the relative un- 
certainty y/V(M)/(M) vanishes in the limit of large N. 
In contrast, there exist "anomalously fluctuating" pure 
quantum states \ip) with a variance V^p(M) ex N 2 . Such 
states are called macroscopic according to Ref. [10J . In 
Ref. [TS], we motivate an assignment of an effective size: 
iV e s ff (^) = maxMUocai V^(M)/N. 

The second measure that we discuss in this context 
is the quantum Fisher information |28| and is related to 



the work of Bjork and Mana. In Ref. [15], we suggest 
to take the usefulness for quantum applications (e.g., in 
terms of quantum metrology) as well as the indication 
of "macroscopic quantum phenomena" (e.g., the possi- 
bility to observe rapid oscillations generated by local 
Hamiltonians) as indications for macroscopic quantum 
states. However, we do not refer to a specific superposi- 
tion as in Eq. (pi) . The usefulness in quantum metrology 
and the possibility of rapid oscillations, respectively, 
can be quantified by the quantum Fisher information 
T (see Ref. [28] and the Appendix and Refs. |15[ [31], 
respectively). Recall that for pure states \(f>) one has 
T{$>i H) — AV^H). Besides the different motivation, 
there hence exists an intimate mathematical connection 
between the Fisher information of pure states and the 
index p, for which reason one can assign the same effec- 
tive size: N^ s (p) = n\ax Hdocal F(p, H) / (AN) (see also 
Ref. H5]). 

Already in the context of the relative Fisher informa- 
tion, we have calculated the maximal variance for local 
operators, V (iJ opt ) = \{N + l) 2 + N + 1, leading to 
an effective size 7V e s ff = iV c F ff = \ N + 2 + 3/(2JV). This 
means that the quantum states producible by means of 
PCQC are macroscopic with respect to these measures, 
see Table |T] We emphasize that the entanglement to a 
microscopic particle is not necessary. In fact, any single 
qubit state that is subject to a PCQC results in a state 
that is in this sense macroscopic. 

Type Authors Effective size 

Bjork and Mana [TT] 1 + 0(1/N) 

MS Korsbakken et al. [TJ] 2 

Marquardt et al. [13] 2 

MQ Shimizu and Miyadera Q2] ±N + O(l) 
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iJV + 0(l) 



TABLE I: Summary of the effective size due to the mea- 
sures of macroscopic superpositions (MS) and macroscopic 
quantum states (MQ), respectively. 

The quantum Fisher information J 7 of the pure cloned 
state \ip) of Eq. {I2J) scales with N 2 , which is optimal for 
local H and implies that the state can be used for quan- 
tum enhanced metrology |32| . However, in the presence 
of local phase noise, the Heisenberg scaling cannot be 
attained anymore, and the best improvement one can 
gain is a constant factor e compared to classical strate- 
gies [33]. While a Greenberger-Horne-Zeilinger state in 
this scenario does not perform better than a product 
state |31j, spin squeezed states [35J can asymptotically 
achieve this limit [35]. Note that Dicke states are not 
spin squeezed [37 . Hence, the simple protocols con- 
sidered in Refs. |34l [35] do not give an improvement 
compared to classical input states. 

Here we study the usability of Dicke states for pa- 
rameter estimation in the presence of different kinds of 
noise. Neglecting the micropart of Eq. pi), we investi- 
gate the performance of the Dicke state Tn, | (N — 1)) 
for a time evolution that consists of a local unitary ro- 
tation around the x axis with frequency w, where in 



addition local bit-flip noise with a decoherence rate 7 
is considered. Note that this is in contrast to Ref. [38 , 
where the use of PCQC was suggested to amplify the 
signal after the unitary time evolution. We consider 
two scenarios: optimal local and global measurements. 
As in Refs. |34|[5rJ] and other references, we give the rel- 
ative improvement of those scenarios compared to the 
best classical protocol. We find that the optimal global 
measurement leads to a better performance than the 
classical strategies, similarly as the spin squeezed states. 
However, within the numerical optimization procedure, 
we were not able to find a symmetric local measurement 
that could beat the classical protocol. This means that 
while full control over the measurement setup directs 
to higher sensitivity, this improvement is lost if only lo- 
cal measurements are available. Qualitatively, the same 
conclusion can be drawn for white noise as a decoher- 
ence source, which is -in addition to illustrations and 
further details- discussed in the Appendix. 



Conclusion and outlook 

We have discussed the generation of macroscopic 
quantum states by means of optimal phase covariant 
cloning. We have argued that neither with proposals 
from the literature nor by arguments on the distin- 
guishability it is justified to call the generated micro- 
macro state (pi a macroscopic quantum superposition 
in the spirit of Schrodinger's cat. It is important to note 
that these results do not depend on the specific map- 
ping of Eq. (fij). They are valid in general for optimal 
phase covariant cloning maps -and hence apply also to 
the experiments of Refs. [5] |T7]- because of the restric- 
tions imposed by the no-cloning theorem. The inability 
of perfect cloning leads to a reduction of the effective 
sizes according to Refs. |11ITT31 [T5] and to a decreased 
distinguishability 

On the other hand, we have seen that this state is 
indeed a macroscopic quantum state since it can show 
nonclassical behavior, in particular an increased sen- 
sitivity in quantum metrology even in the presence of 
noise. Our results indicate that a careful investigation 
is required to judge whether large-scale quantum sys- 
tems can indeed be considered to be truly macroscopic. 
While we have concentrated here on systems consisting 
of many spins, similar investigations, e.g., for superpo- 
sition states produced in optomechanical systems are 
highly interesting, given the spectacular experimental 
progress in this area. 
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Measuring a subgroup 

To determine the effective size of the macro-measure 
of Korsbakken and coworkers [IS], we calculate the dis- 



tinguishability of \ip + ) and \ip~) if we only measure k 
out of N particles (see also Ref. |B|). The figure of merit 
is the success probability P = | +■ \\\p^ — p% ||i with 
p^ = Tr N \ k |i/' ± )(V' ± l- We use the general decomposi- 
tion 
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where ( ) is the binomial coefficient that gives the 
number of possible permutations. Defining dj — 
we get 
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We consider in the following large particle numbers N. 
For x € {^(^ — 1), |(JV + 1)} we can approximate 
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( N x Z-)/0 = 2-*[l + 0(*/JV)]. With ^c—cj^ 
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c 2 y/(k — j)/(j + 1), the trace norm of 
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(\k,j){k,j + l\+h.c.) 



can be easily computed numerically. From Eq. (HI, the 
same calculations can be done to find p N _ k in the limit 
of large N . The results are discussed in the main text. 



culation of C renders easily and we get 
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The matrix consists of three independent 2 x 2-blocks: 
an identity matrix and two times the matrix 
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JV+1 
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JV+1 
JV 2 +2JV-1 



which belong to the coefficient pairs a x 
and {(^i,y,OL2,y) 1 respectively, 
the maximal expectation value 



We 
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— (Xi,x<X2,x {N + !)■ After inserting 



the normalization condition |ai )X | = y/N (1 — ct2, x ) + 1 
(and fixing a\ tX to be positive), we find that a x = 
(1,-1) is an optimal choice. This leads to a maxi- 
mal variance V^,(Af opt ) = |(iV + l) 2 + N + 1. No- 
tice that any collective local rotation around the pr- 
axis e (-^z)®(N+l) M opt e (ia<r z )®(N+l) ^ a g R giyes the 

same variance. 

In addition, we search for the optimal measurement 
for the |</>o) = |+) <8> |"0 _ ) state. Since it is a prod- 
uct state between the first particle and the rest, the 
variances of both parts sum up V,p (M) — V\+\(M) + 
Vw>_ (M) . The micro-part has with a y a maximal vari- 
ance equaling one, while for the macro-part one can use 
once again the covariance method for the measurement 
M = ^2j =x y z Oj,2 X,;=2 a j ■ This leads to the matrix 



N 2 



Maximal variance of micro-macro state 



Given the N + 1 particle state \<f>), we search for 

the local operator M = J2i=i Ei=i y z a i.j a j w ^h 
maximal variance V<j,(M), under the condition that 

EiLi'E™,*^ = N + 1. We need this in or- 
der to decide on the effective size of several mea- 
sures [TtJl [TT1 [T5] . The covariance-matrix method [2T] 
is an efficient method to find the optimal set {otij}. 
With Act:- = <7 — (cr- ),/, we define the hermitian 
covariance-matrix [(«,j) now serve as a multi-index 






(AcrjAaj)}^,]. The maximal eigenvalue 
of C is the maximal variance of a local observable M 
and entries of the corresponding eigenvector are -up to 
a normalization- the weights ctij. 

For \ip) in Eq. pi of the main text, the dimen- 
sion of C is even more reduced, since the macro-part 
of \ijj) with N particles is symmetric under particle 
permutations and therefore we can make the ansatz 



M = T,j=x, y ,z a JA< T j +«J,2E i= 2 °i • The normal- 
ization reads Y] . 



y,z a JA 



a 2 2 N = N + 1. The cal- 
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1, 



and shows that the maximal variance is governed 
by a measurement in y-direction, which leads to 
V^M^) = 1 + \{N 2 +2N - 1) = ±(JV + l) 2 . 



Details on the distinguishability 

This paragraph contains some details of the calcu- 
lation concerning the distinguishability of \ip + ) and 
\ip~) under imperfect conditions. We therefore calcu- 
late the probabilities (7Tj)^± which are used in the first 
scenario in the main text. Using Eq. (El, the states 
|-0 ± ) can be expressed as superposition of Dicke states 
\N,k) x = Had® \N,k) in the :r-basis (Had denoting 
the Hadamard operator) 
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N,k) x (9) 



with a, = 1/^-1(40 Et (-ir©(4i-,)- 

With this we easily see that 
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(10) 



In the same manner, one can calculate the proba- 
bilities under the assumption of a POVM with Kraus- 
operators Ej defined in the main text. We arrive at 



iV 



(Ef 



v [i±(-i)T 

^ 2 



fc=0 



Ptn k e~ 



(11) 



A numerical example for (Ef)^,± is shown in Fig. pi 

Exchanging the pure input states l^) by noisy 
£z (IV' )(V' ± I) is more involved. An error on a spin par- 
ticle leads to a wrong assignment of the measurement 
outcome, similarly as with the POVM measurement, 
but multiple errors can cancel this effect. Considering 
carefully all posibilities, we end up with the formula 



pf =Ik[7r i £: i (|V^XV^|)] 

» \i±(-iyf (N\ 

= 1^ 2 I «■ > 0kVk 

fe=0 ^ 



(12) 



with 



N/2 
, li-fel 



.N-21 



"-«^ l+ % 



N-k\ 

;_*_-,)• (13) 



Fig. [3] nicely illustrates the discussion in the main 
text. The distinction can be done on two different 
scales: on the "microscopic" scale, where we can per- 
fectly distinguish between \tp ) measuring the parity, 
and on the "macroscopic" scale, where one tries to com- 
pare the total number of excitations. The micro-level 
is very fragile against any kind of disturbance and the 
distinguishability is lost almost instantaneously. The 
more robust distinction on the macro-level suffers from 
the large variances of the observable and one therefore 
cannot distinguish with arbitrary precision. 

Parameter estimation protocol 

Here we give some details on the time evolution and 
the measurements to calculate the so-called Cramer- 
Rao bound [5S] for the quantum state \N, |(iV — 1)). 
The Cramer-Rao bound is a lower bound on the mini- 
mal error Su> one has to assign for the estimation of a 
parameter u. The time evolution that we consider is 
governed by the master equation 



m 



N 



*£<# 



i=l 



N 



(=1 



(14) 



with the parameters u > for the unitary rota- 
tion and 7 > for the decoherence rate. The 




FIG. 3: Probabilities p, - {Ef)^± are plotted for N - 31 
and thre e different choices of a: [i.e., perfect measurement 
Eq. 1 10 1], 1 and y/N (from top to bottom). Already a uncer- 
tainty of a — 1 reduces the distinguishability between |V-' + ) 
(blue) and \4>~) (red) drastically. 



quantum state under this evolution is pit) = 

e (-tort/2)eJV£ B (1^ N-ly^ N ^N-l^ e ( 4wt /2)8JV ] recalling 

that S x {p) — Yli=i U P + (1 — M ) cr xpcrx with u — 
i(l + e-^)e[|,l]. 

We calculate the Fisher information of p(t) for the 
optimal measurement -the so-called quantum Fisher in- 
formation J-- and a specific local measurement Mj = 



N 



,w 



The latter is also referred to as classical 



Fisher information F. Here, T is given by 



F(p,H) = J2 2 



Pi(t)+Pj(t) 



mtH/2\j)\ 2 , (15) 



where Pj(t) are the eigenvalues of the density matrix 
p(t), that is, the evolving state without unitary rota- 
tion. The corresponding eigenvectors \j) are time in- 
dependent. We calculate the spectral decomposition 
of p(t)\ u =o numerically and are therefore restricted to 
rather small system sizes. 

The the Cramer-Rao bound in this scenario reads 
5uj = 1/ ' \fnJ- ' . where n is the total number of repe- 
titions. Assuming a total time T for the experiment, 
we have n = T/t. We then minimize 5u with re- 
spect to t. We compare dui to product states as in- 
put states. The same procedure leads there to the un- 
certainty Sujps = y/2e r y/(TN) [33] • We are interested 
in the relative improvement 1 — Suj/Sujps, which is the 
quantity plotted in Figs. [4] and [5j 

Similarly, we proceed with the classical Fisher in- 
formation F. For a discrete probability distribution 
Si(uj,t) that depends on the parameter u, F is defined 



by 



* = £• 



d log a, 

d'jj 



(16) 



Here, s,(w, t) are probabilities of the results for the dif- 
ferent outcomes of a given measurement M,. We start 



with the measurement M, 



= Ei=iO" 



W 



where we can 



find analytical expressions for Si(oj,t) 



I pit) 



where 



7Tj is the projector on the subspace spanned by all prod- 
uct states with i times |1) states and N — i times |0) 
states. This is done by using the splitting in Eq. (EJ for 
|JV, \{N - 1))(N,\(N - 1)| with respect to i : N-iand 
a consequent expression of the sub-parts in the x-basis, 
similar to Eq. @. We find 



Si(u,t) = (ni) p (t) 
= ( N \ (1®* ® O^-'l p(t) |l®* <8) O^' 4 ) 

E-pi.- r N-i.+ 

3d'=J 



(17) 



with the definitions 



A*j = 



(")Q)(4- r -,) 

N (40 



r S = E (±i) w H*4,^ 



i,i/,m=0 

z— m 7?i —iu)t(l-ll) 
U U e v ', 



'J,l 






£(-*) 



m=M 



ml "mil' 



i-l 

j -m 



(18) 



In Eq. (171 


we have J = 


{max[0,i- |(JV + 1)) , 


. . , min ( | (JV + 1) , i) } and in 


Eq. (181 M = {max(0,j 


-i + /),..., min(Z,j)}. 



From this we can calculate F and pursue as before 
with J- '. To improve the results, we can try to opti- 
mize the local measurement. We test local measure- 
ments where every qubit is measured in the same basis. 
Furthermore, we assume to find the best basis in the 
y — z-plane of the Bloch sphere of a single spin, because 
measuring in x-direction only gives time-independent 
results. Instead of |0) and |1), we now project onto 



|0) and e 



|1), a e [0,7r/2]. Thi s can be 



easily implemented in Eqs. (17 1 and (18 1, exchanging 



the term e ~^Ki-i') by e -i(^Wii){l-U) ^ We 8earc h nu 
merically for the optimal a and present the results in 
Fig. [4] We find indeed an improvement with respect to 
M z , which is however still above the uncertainty <5wps- 
Whether a more general, asymmetric measurement can 
improve the sensitivity is an open question. 

The very same can be exercised with different noise 
channels, as long as it commutes with the unitary 
evolution in a;-direction. An example is the so-called 
white noise. This channel is defined as £w(p) = 

The 



HL{pp- 



^ J2i=o °iPOi with p = 



= p-rt 



results for the quantum Fisher information ( 15 1 are 



slightly better compared to the phase-error channel £ x 



while the formulas Eqs. (17 1 and (18 1 are indeed also 
valid for white noise if we identify u = (1 + p)/2 (see 
Fig. [5). 
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FIG. 4: Relative improvement of phase estimation proto- 
cols using \N, I (N — 1)) as a starting state with the param- 
eters w = 1 and 7 = 0.5. We compare the three scenarios 
global and local optimal measurement as well as fixed mea- 
surement in z-direction to the best strategy with product 
states as inputs. The global measurement can lead to an 
increased sensitivity for the investigated system sizes. Note 
that the maximal possible improvement is 40% |33| . On the 
other hand, local measurement strategies fail to overcome 
classical protocols. The different ranges of data points re- 
sults from different numerical techniques (see text). 



■ Optimal global measurment 

■ Optimized local measurement 
. M -measurment 




15 20 

System size N 



FIG. 5: With white noise as decoherence source, the uncer- 
tainty of the estimation protocol compared to classical strat- 
egy is the same for our choice of local measurements, while 
the relative improvement changes for the optimal global 
measurement (see also Fig. B; parameters: u) — 1 and 
7 = 0.2. 
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